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Numerical evidence of the axial magnetic effect 
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The axial magnetic field, which couples to left- and right-handed fermions with opposite signs, 
may generate an equilibrium dissipationless energy flow of fermions in the direction of the field even 
in the presence of interactions. We report on numerical observation of this Axial Magnetic Effect in 
quenched SU{2) lattice gauge theory. We find that in the deconfinement (plasma) phase the energy 
flow grows linearly with the increase of the strength of the axial magnetic field. In the confinement 
(hadron) phase the Axial Magnetic Effect is absent. Our study indirectly conflrms the existence 
of the Chiral Vortical Effect since both these effects have the same physical origin related to the 
presence of the gravitational anomaly. 
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Anomalies belong to the most characteristic and fun- 
damental properties of relativistic quantum field theories. 
They signal an incompatibility between quantization and 
the symmetries present at the classical level. While the 
effects of anomalies in vacuum are well understood it has 
only recently been fully appreciated that anomalies play 
also an extraordinary important role at finite temper- 
ature and density. In particular they give rise to new 
non-dissipative transport phenomena. The most well- 
known of these is the so-called Chiral Magnetic Effect 
(CME) [1], describing the generation of an electric cur- 
rent parallel to a magnetic field in the presence of an 
imbalance between the number of right-handed and left- 
handed fermions. The CME is thought to be responsible 
for charge asymmetries observed in heavy ion collisions at 
RHIC and LHC [2 . It also might play a role in the trans- 
port properties of advanced new materials, the so-called 
Weyl semi-metals in which the effective charge carriers 
can be modeled as 3 -f 1 dimensional Dirac fermions [5] . 

The CME is however only one representative of a whole 
class of anomaly related transport phenomena. A full 
classification of such phenomena has been obtained via 
Kubo formulas in f?|. It turned out that not only the 
usual axial or chiral anomalies give rise to dissipationless 
transport but that there is also a distinguished place for 
the axial gravitational anomaly. 

In general, anomaly related transport is sourced by ei- 
ther external magnetic fields or by vortices in the fluid 
of chiral fermions [5rj Thus we can distinguish between 
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An early manifestation ot this ettect in rotating ensembles of 



chiral magnetic and chiral vortical effects. The gravita- 
tional anomaly comes in through the chiral vortical effect 
(CVE). Even in the absence of chemical potentials the 
gravitational anomaly gives rise to a chiral vortical effect 
at finite temperature 

Js^crw, (1) 

where J5 is the axial current and a; = V x t; is the vorticity 
of the fiuid velocity v. In the absence of matter, the 
conductivity 

depends on the temperature T and the gravitational 
anomaly coefficient. Equation ([T]) is valid for a theory 
consisting of massless fermions. In a basis of left- and 
right-handed Weyl fermions qi (g^) are the charges of 
the left-handed (right-handed) fermions. This effect has 
been confirmed at strong coupling via the gauge-gravity 
correspondence in Ref. . 

The chiral vortical effect ([T]) is quite remarkable since, 
contrary to the CME, it does not rely on the explicit 
introduction of an axial chemical potential: in Eq. ([T]) 
the difference of the right- and left-handed charges in 
Eq. ([2]) corresponds to the vacuum content of the theory. 
On the other hand, its relation to the axial gravitational 
anomaly is at the moment still somewhat less direct than 
the relation of the CME to the axial anomaly. Purely hy- 
drodynamic arguments such as in Ref. [H] are not able to 
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fix the prefactor in the conductivity Q due to a mis- 
match in the order of derivatives in which the gravita- 
tional anomaly can apparently influence hydrodynamics 
(see however [ST). 

It has further been proven that the transport law ([T]) 
does not get renormalized in perturbation theory in theo- 
ries which contain fermions and scalars only. If there are 
however dynamical gauge fields that contribute to the ax- 
ial anomaly then it was shown that a non-vanishing two 
loop contribution arises [10]. In QCD this is of course 
the case. The axial anomaly has a gluonic contribution 
and therefore one expects a strong renormalization of the 
chiral vortical effect in QCD. 

Since the CVE is a dissipationless and stationary ef- 
fect it is accessible via Euclidean field theory and can, 
in principle, be studied numerically in simulations of lat- 
tice gauge theories. A straightforward implementation of 
rotating fluid on a lattice seems, however, a rather non- 
trivial task due to obvious incompatibility of the small, 
discrete rotational symmetry group of the Euclidean lat- 
tice with smooth, continuous rotations used in Eq. ([T]). 
Fortunately there is an alternative way of accessing this 
particular transport phenomenon that is well suited for 
implementation on the lattice. 

In order to compute the CVE for the axial current one 
could make use of the Kubo formula 

where J5 are spatial components of the axial current 
and T^'^ are temporal-spatial components of the energy- 
momentum tensor. Since the correlator is to be evaluated 
at zero frequency one can reverse the order of the opera- 
tors J| and T^*^ in Eq. ([s]) and obtain a new effect corre- 
sponding to the generation of an energy current = T'^^ 
in the background field that couples to the axial current, 
this is an axial magnetic field. One finds therefore that 
the chiral vortical conductivity ^ also appears in the 
new transport formula ^ : 

I = 'yB^ , (4) 

which represents the Axial Magnetic Effect (AME). The 
transport law Q describes the generation of an equilib- 
rium dissipationless energy current in the presence of an 
axial magnetic field at finite temperature. 

The practical advantage of the AME formula Q is 
that the axial magnetic field i?5 can be relatively easy 
implemented on the Euclidean lattice, while the imple- 
mentation of the vorticity oJ is a much more difficult 
task. On the other hand, both the CVE and AME have 
the same physical nature - which is also clear from the 
very fact that they share the same conductivity coeffi- 
cient ^ - originating due to the presence of the gravita- 
tional anomaly [4 . Thus, in this article, we concentrate 
on numerical evaluation of the AME law ^ in the con- 
text of the quenched SU(2) lattice gauge theory for three 



different temperatures, which represent three basic re- 
gions of the phase diagram: the deconfinement regime, 
the critical confinement-deconfinement region, and the 
confinement phase. 

The important feature of the AME is that it is realized 
in the pure vacuum with all chemical potentials equal to 
zero, fi = fi^ = 0. The dissipationless equilibrium energy 
fiow is achieved at finite temperature, T 7^ 0, in the pres- 
ence of the axial magnetic field B5, which distinguishes 
left-handed and right-handed quarks. The axial magnetic 
field couples to the left-handed and right-handed quarks 
with opposite charges respectively. 

In order to check the existence of the AME law Q it 
is sufficient to consider one type of fermion with a unit 
charge: 

q5=-'li^+e- (5) 

The coupling of the quarks to the chiral field is de- 
scribed by the following Lagrangian: 

£5 = ^{d^ - tgA-^t- ~ ij'eA5,^)r^P = 4>'PM^)^ , (6) 

where is the nonabelian SU (N) gauge field and are 
the generators of the SU{N) gauge group. 

The energy flow in the transport law Q is given by 
the expectation value of the off-diagonal component of 
the stress-energy tensor T'^'^, 

Jl^iT'n^lii'il^'Dl+Y'DlW, (7) 

where the latin index i = 1, 2, 3 labels the spatial coordi- 
nates and = is the time direction. 

Wc introduce the stationary uniform axial magnetic 
field in the third direction B^^i — ■ Si^s by setting 
^5,0 = ^5,0 = 0, ^5,1 = -X2B5/2 and ^5^2 = X1B5/2. 

The energy flow can be implemented on the lattice via 
a straightforward discretization of Eq. ([t]) using the linear 
combinations of the nonlocal lattice correlator, 

C^{x,y;A) = {^(x)U,^y{A';,h^i;{y))A 

In this formula, the expectation value is taken over the 
fermion fleld in a flxed background of non-Abelian 
and axial A^^^ flelds, and the trace is taken over color 
and spinor indices. Here Ux,y is the gluon string between 
the lattice points x and y which makes Eq. ([S]) gauge 
invariant. The Dirac operator I?5 is given in Eq. (§). 

We calculate the correlation functions ([s]) numeri- 
cally, using lattice Monte-Carlo simulations of quenched 
SU{2) lattice gauge theory following numerical setup of 
Refs. [TT]. The quark fields are introduced by the overlap 
lattice Dirac operator V with exact chiral symmetry [T^] . 
The correlation functions ([s]) are substituted in the dis- 
cretized version of Eq. Q and then the whole expression 
is averaged over an equilibrium ensemble of finite tem- 
perature configurations of non-Abelian gauge fields 

(O) = (^j DAI j DAI e"^™^^"'! O , 
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where S'ym(^^) is the lattice action for the gluons 

There are at least two ways to calculate the fermion 
propagator (|8| in the presence of the axial magnetic field 
in a finite-volume lattice. One can introduce the axial 
magnetic field straightforwardly by modifying the spatial 
boundary conditions for fermions according to the gen- 
eral approach of Ref. [17. Alternatively, one can make 
use of the identity, 



tr [55(^5) 7 J = tr [{Vr + VL)SdA5) 7^ 
= tr [Vr 5(^5)7,,] + tr [Vl Si-A^)^^] , 



(9) 



where Vr,l = (1 ± 7^)/2 are the right and left chiral 
projectors, the trace is taken over spinor indices and 

55(^5) = [1^5(^5)] , S{A) - [TpiA)] , (10) 

are the Dirac operators for the massless fermions in the 
background of the axial field vls.^ and the usual Abelian 
gauge field A, respectively. The former is defined in 
Eq. (|6| while the later has the usual form: 

Vf,{A) = df,~igAp''-ieA^. (11) 

It worth stressing that in the right hand side of Eq. ([9| 
the axial gauge field ^5 appears as the Abelian field with 
opposite signs for right-handed and left-handed fermions. 
This is an expected property given the very definition of 
the axial magnetic field. 

Identity ([9| is valid regardless of the dynamical gener- 
ation of quark mass and chiral symmetry breaking since 
this identity is a generic property of the fermion operator 
itself while the mentioned phenomena are the particular 
properties of the expectation values of this operator. 

Thus, identity ([9| allows us to express the energy 
fiow ([7]) of the massless fermions via the standard Dirac 
operator in a background of usual magnetic field. This 
property is particularly useful for numerical calculations 
with overlap fermions since we can use the already exist- 
ing techniques of Ref. [11] . 

We evaluate the energy flow ([?]) in the quenched SU (2) 
gauge theory using 300 configurations of the gluon gauge 
field for each value of the background axial magnetic field. 
We consider the asymmetric lattices L^Lt with three 
temporal lengths Lt = 4, 6, 8 and the fixed spatial length 
Ls = 14. We use the improved lattice action for the 
gluon fields with the lattice coupling (3 — 3.2810 which 
corresponds to the lattice spacing a — 0.103 fm [14]. 

Similarly to the usual magnetic field, the axial mag- 
netic field is quantized due to the periodicity of the gauge 
fields in a finite volume : 



S5 ^fcSs, 



2n 
if 



0.117 GeV^ (12) 



where the integer A: = 0,l,...,-Lg/2 determines the num- 
ber of elementary magnetic fluxes which pass through the 
boundary of the lattice in the {x^,x'^) plane. Notice that 
the elementary (minimal) field ( 12 ) in our simulation is 
three times weaker compared to the field used in our pre- 
vious studies, Ref. [llj . because in the present paper the 
fermion field carries a unit electric charge. 



The maximal possible value of the quantized flux, k = 
L^/2 = 98, corresponds to an extremely large magnetic 
field. In this case the magnetic length Lb ^ (eS)^^/^ 
is the order of the lattice spacing, Lb ^ a. In order to 
avoid possible ultraviolet artifacts, we consider relatively 
weak axial magnetic fields with k = 0, 1, . . . fcmax, where 
the fiux number is limited by fc,„ax = 10 <C ^^/2, so that 
our strongest magnetic field is eiJmax ~ 1.17 GeV^. 

In order to increase the efficiency of our numerical algo- 
rithm we introduce a small bare quark mass m ^ 20 MeV 
in the overlap fermion operator Q. Despite this bare 
mass being very small, we have carefully checked the ap- 
plicability of Eq. ^ to our numerical setup by making 
sure that the energy fiow is insensitive to the variations 
of the quark mass. For example, a two-fold increase of 
the mass leads to a less than 1% change of the central 
values of our observable (this is much smaller than our 
statistical errors). 

In SU{2) gauge theory the critical deconfinement tem- 
perature is Tc = 303 MeV. Our three lattices with tem- 
poral extensions Lt — 4, 6, 8 correspond, respectively, 
to the deconfinement phase (T = 1.58 T^), the vicin- 
ity of the confinement-deconfinement phase transition 
(T ^ 1.05 Tc) and the confinement phase (T = 0.79 T^). 
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FIG. 1. Energy flow ^ parallel (j|) and perpendicular (Je*") 
to the direction of the axial magnetic field £5 in the decon- 
finement phase at T = l.SSTc = 479 MeV. The red dashed 
line represents the best linear fit confirming the existence of 
the Axial Magnetic Efi^ect (j4|. The slope of the fit is given by 
the conductivity (13 1. 



In Fig. [T] we show the energy flow ([t]) as a function of 
the axial magnetic field in the deconfinement phase 
at T = 1.58 Tc. One can see that the energy flow parallel 
to the axial magnetic fleld j]' is a linear function of the 
field strength as predicted by the AME transport law Q . 
The energy flow in the transverse directions is zero 
within the error bars, as expected. We have also checked 
numerically that the usual magnetic fleld does not induce 
any energy flow. 

The linear fit of the energy flow gives us the following 
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value of the conductivity: 



aiT) 



T=1.58Tc 



= 2.22(3) X 10"^ GcV^ 



(13) 



magnetic field. Our study also indirectly confirms the 
existence of the Chiral Vortical Effect (CVE) since both 
these effects have the same physical nature originating 
from the gravitational anomaly. 



For a conformal theory, the conductivity coefficient is 
expected to be crthcor — T^/12 which gives us one order 
of magnitude larger value trthoor ~ 1.9 x 10-2 GeV^ at this 
temperature. This difference is not unexpected because 
at T = 1.58 Tc the Yang-Mills theory is stiU far from the 
conformal regime of high temperatures. 

We see no signature of energy flow neither in the vicin- 
ity of the phase transition nor in the confinement phase. 

In the high temperature limit the theory is expected to 
become conformal and the conductivity a{T) should be 
described by Eq. Given our results, we suggest that 
in QCD the temperature behavior of the conductivity 
may be described by Eq. ^ multiplied by a prefactor 
/(T) which approaches unity in the high temperature 
limit while damping (likely, exponentially) the conduc- 
tivity as the temperature gets decreased. 

Summarizing, we have numerically observed the ex- 
istence of the Axial Magnetic Effect (AME) in the 
quenched SU{2) gauge theory on the lattice. We 
have found that the equilibrium energy fiow of massless 
fermions is parallel to the direction of (and proportional 
to the strength of) the time- independent uniform axial 
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